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ABSTRACT 

We present phase-resolved optical spectroscopy of the counterpart to the neutron star 
low mass X-ray binary EXO 0748—676, almost one year after it turned into quiescence. 
The spectra display prominent H/3 and H7 and weak Fe 11 lines in emission. An average 
of all the spectra (corrected for the orbital motion) also exhibits a very weak line from 
Mg I. Tomographic reconstructions show that the accretion disc is not contributing to 
the optical line emission, which is instead dominated by the irradiated hemisphere of 
the companion star facing the neutron star. We could not detect absorption features 
from the mass donor star in the spectra. The emission lines appear broad, with an 
intrinsic FWHM of 255±22 kms~^. Under the assumption that the width of the 
Fe II emission lines is dominated by rotational broadening, we obtain a lower limit 
on the compact object mass which is inconsistent with a NS accretor. We discuss 
this incongruity and conclude that either the lines are blends of unresolved features 
(although this requires some fine tuning) or they are broadened by additional effects 
such as bulk gas motion in an outflow. The fact that the Fe 11 lines slightly lag in phase 
with respect to the companion star can be understood as outflowing gas consistent 
with a Black- Widow like scenario. Nevertheless, we can not rule out the possibility 
that blends of various emission lines cause the apparent phase lag of the Fe il emission 
lines as well as their large width. 

Key words: stars: individual (EXO 0748—676) — accretion: accretion discs — stars: 
binaries eclipsing — X-rays: binaries 



1 INTRODUCTION 

EXO 0748—676 was discovered with the European X-Ray 
Observatory Satellite (EXOSAT, |Parmar et al.||1985| ) in 
1985. Soon after the discovery, the detection of type I X- 



ray bursts from the source (Gottwald et al. 19861 marked 



it as a Galactic low-mass X-ray binary (LMXB) where a 
neutron star (NS) is accreting matter from a low-mass com- 
panion star. Many LMXBs are known to be transient, al- 
ternating periods of quiescence at a relatively low X-ray lu- 
minosity (~10'^'^ ergs~^in a 0.5-10 keV range) with month 
to year-long bright X-ray outbursts (10'*'' — 10^* ergs^^). 
EXO 0748—676 has been continuously in outburst for the 
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24 years since its discovery: tens of years-long X-ray out- 
bursts have been observed from a number of LMXBs, e.g. 



KS 1731-260 ( |Wijnands et al.||2001[), G RS 1915-M05 and 
4U 1755-338 (sec Rcmillar d fc McCh ntock 20061 for » re- 
view). During outbursts mass accretion proceeds at a high 
rate forming an extended accretion disc around the accreting 
compact object. X-rays are emitted from the innermost re- 
gions, whereas the optical originates further out in the disk. 
In LMXBs the disc dominates the optical flux during out- 
bursts, outshining the low-mass companion star. The latter 
can become visible during quiescence, when the disc is less 
bright. 

A few LMXBs are known where the NS is detected as a 
pulsar ( |Chakrabarty|2005| but for the majority of them the 
companion star is the only viable source of information re- 
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garding the system dynamics. Optical spectroscopy of the 
companion star can be used to measure the orbital param- 
eters and, under certain conditions (see below), the mass of 
the accreting compact object. This measure is interesting 
in LMXBs since the accretion process can significantly in- 
crease the mass of a NS: the maximum mass that a NS can 
reach is one of the parameters that can distinguish among 
the existing models for the equation of state (EoS) of those 



object (Lattimer & Prakash 20011. The determination of 
the NS EoS is one of the key goals in the study of NSs and 
will have strong implication for both astronomy and super- 
nuclear density matter physics. 

The mass of a NS (or a BH) in a non-pulsating LMXB can 
be measured by solving the system mass function 
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where G is the universal constant of gravity, A'Id and Mx 
are the masses of the companion and NS, respectively, P is 
the orbital period, i the inclination of the orbital plane, Kd 
the amplitude of the radial velocity curve of the companion 
and where we define q = Md/Mx- This equation is valid 
provided that Kd represents the motion of the center of 
mass of the companion. Because /{Md) ^ Mx, the mass 
function provides a lower limit to the mass of the neutron 



star ( Charles & Coe 2006 1 



P and K_D can be inferred from the orbital Doppler shift of 
stellar absorption lines originating in the atmosphere of the 
companion star, which can be visible in the optical spectra 
during quiescence. /(M_d) can be solved for Mx in eclipsing 
quiescent LMXBs ( |van Paradijs fc McClintock|1995l ), where 
q can be expressed as a function of Kd and of the projected 
rotational velocity of the companion star v sin i = ^Rd sin i 
| |Wade fc Hornej 198 81 measured from the broadening of the 
companion stellar absorption lines ( Gray|1992_) . From q and 
from the eclipse duration, i can be obtained |^ 



The LMXB EXO 0748-676 shows X-ray ecUpses (Parmar 



et al. 1986 1 and, after a more than 20 year-long outburst, 



it turned into quiescence in September 2008 (Wolff et al. 



2008b| [Wolff et al.|2008al|Hynes fc Jones|2009, Torres et al. 



2008). The optical counterpart to the X-ray source was first 



found by Parmar et al. ( 1985 \ and confirmed by the Chan- 
dra localisation of E XO 0748-676 at RA= 07''48'"33':73, 
Dec.= -67°45'07'.'9 (iTorres et al.||2008l). A search on pho- 



tographic plates showed that the source was not detected 
down to ~23 magnitudes when the X-ray emission was off. 
Spectroscopic observations in the optical during the out- 



burst have been performed by Pearson et al. ( 2006 1 and 



Munoz-Darias et al. (20091, attempting to get a lower limit 



on the NS mass from He, C and N emission lines probably 
originating from the inner, heated face of the companion 
star. The first spectroscopic study of EXO 0748—676 in qui- 
escence was performed by Bassa et al.| ( |2009[ ) two months 
after the end of the outburst. The spectra, acquired in a 
5750-7310 A wavelength range, were dominated by strong 



sin-^ i cos^ (tt At/)) = 1 — [- 



I 499^3 ^ 

eclipse duration (see |Horne|1985| . This relation is valid under the 
assumption that the companion is Roche Lobe filling, which holds 
for LMXBs. 




0.6 

Orbital phase 

Figure 1. g'-band light curve determined from the acquisition 
images. The magnitudes are relative to the bright star USNO- 
Bl.O 0222-0189796. Different symbols refer to different observing 
nights starting from the night of 2010 January 18/19. In tempo- 
ral order, the observations of the first to the fourth nights are 
indicated by black stars, black squares, empty circles and empty 
squares respectively. 



Ha and weaker He emission lines coming from the compan- 
ion star facing the NS, an indication of irradiation. A weak 
contribution to the lines was due to optical radiation from 
a residual accretion disk. No absorption lines were detected, 
preventing a straightforward solution of the f{Mo)- Never- 
theless, the authors could put a lower limit on Kd and on 
the NS mass, Mns > 1.27M0. 



Bassa et al. ( 2009 1 determined a minimum temperature for 
the irradiated companion star of EXO 0748-676 of ~5000 
K, the surface temperature of a G type star. We have per- 
formed phase- resolved spectroscopic observations of the op- 
tical counterpart to EXO 0748—676 after almost one year 
of quiescence. We attempted to detect absorption lines from 
the heated companion star by collecting our spectra in the 
wavelength range (4222-5701 A) where a G as well as a A" 
type star would show absorption lines. We here report on 
the results of this study. 



2 OBSERVATIONS AND DATA REDUCTION 

We performed long-slit phase-resolved spectroscopy of the 
optical counterpart to EXO 0748-676 with the FORS in- 
strument on the Very Large Telescope (VLT) (grism 
1200g-(-96 and a I'.'O slit). The orbit was sampled with a 
total of seventythree 900s long exposures collected between 
2010 January 18 (MJD(UTC)=55214) and 2010 January 22, 
in the wavelength range 4222-5701 A. We also acquired tem- 
plate spectra from main sequence G type stars (G5V, G9V 
and G6V). The seeing across the four nights was between 
~ 0'.'8 and ~ l'.'2. The detectors were read out with a 2x2 
binning, providing a resolution of 2 A (measured from the 
width of both arc lines and of the night sky OI line at 
5577.338 A) sampled with a dispersion of 0.73 A pix"^ The 



VLT observing program 085.D-0441(C) 
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Figure 2. Trailed spectra composed by 73 single spectra of the counterpart to EXO 0748—676, phase binned in 20 bins. The two 
prominent s-waves are from the H7 and H/3 lines. Three weak lines consistent with Fe II are just visible between 4900 and 5200 A (see 
section |3.1| l. 



images were corrected for bias, flat-fielded and extracted us- 
ing the Figaro package within the Starlink software and 
the packages Pamela and Molly developed by T. Marsh. 
We used dome flats for the flat-fielding and we subtracted 
the sky continuum by fitting clean sky regions along the 
slit with a second order polynomial. The spectra were opti- 
mally extracted following the algorithm of [Home] ( |1986| im- 
plemented in Pamela and wavelength-calibrated in Molly 
with a final accuracy of 0.1 A, using arc exposures taken 
during daytime. The wavelength calibration was corrected 
for shifts in the single observation with respect to the po- 
sition of the sky OI line at 5577.338 A (Osterbrock et al 



1996). Each spectrum has been normalised dividing by a 



first-order polynomial fit of the continuum. The spectra have 
been phase binned with the Tp = 54776.501663 ± 0.000068 



MJD/TDB from Bassa et al. (20091, which is the closest 



in time to our observations and with the orbital period 
Port = 0.15933783446 days fromlWolff et al.l ([20091. 



We also analysed sixteen 6 s-long g'-band acquisition im- 
ages, corrected for bias and fiat-fielded with standard rou- 
tines running in MIDAS. The photometry was performed 
through point spread function fitting, using daophot II 



(Stetson 19871. Absolute photometric calibration was not 



possible due to the lack of observations of g'-band stan- 
dard stars, but a light curve (Fig. [T]) of the optical coun- 
terpart to EXO 0748—676 was obtained through relative 
photometry with respect to the reference star USNO-Bl.O 
0222—0189796. The instrumental magnitude of the reference 
star was measured with an accuracy of 0.03-0.08 magni- 
tudes across the various images. Relative photometry with 
respect to two other bright targets have shown that the ref- 
erence star was not variable during our observations. The 
light curve was phase-folded with the same ephemeris used 
to phase-bin the spectra (see above). 



3 ANALYSIS AND RESULTS 

3.1 Individual and trailed spectra 

The sole features detected in the individual spectra are the 
Balmer lines H/3 (4861.327 A) and H7 (4340.465 A), which 
are observed as strong emission features. We searched for ab- 
sorption lines in the individual spectra by cross-correlating 
with template star spectra, but we found no correlation. 
As shown in Figure [2] more emission features become ob- 
servable when trailed spectra are constructed. Besides the 
s-waves associated with the Balmer lines, three weak s-waves 
are visible in the region between 4900-5300 A: the position 
of the first two lines is consistent with a couple of He I lines, 
at 4921.929 A and 5015.675 A, but also with Fe II hues at 
4923.92 A and 5018.44 A. No Hel line matches the posi- 
tio n of the third line, which instead can be Fell at 5169.03 
A ( [Moore[[l972 |. This favors an interpretation of the three 
lines as an Fe triplet as all the three lines are part of multi- 
plet 42. 

The variations show only one peak along the orbit. The emis- 
sion lines are fading away and disappearing between phase 
~0.75 and ~0.25, consistent with an emission region located 
on the inner face of the companion star. There is no indica- 
tion of absorption lines in the trailed spectra. 



3.2 Radial velocity of the Balmer lines 

For the H/3 and H7 lines we have obtained radial velocity 
curves (see Figure |3| by fitting single Gaussians to those 
lines in each 900 s spectrum, including only significant de- 
tections on a 3cr level. The fitted parameters were normal- 
isation, full-width-at-half-maximum (FWHM) and velocity 
offset of the centroid with respect to the rest frame line 
wavelength. The 0.1 A la uncertainty on the velocity off- 
sets due to the wavelength calibration has been added in 
quadrature to the error on the centroid from the Gaussian 
fit of each observation. We fitted the radial velocity curve 
(velocity offset versus phase) for each line with a circular 
orbit in the form v(<^)=7-|-KcmSiw(27r(/> -j- ip). We measured 
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7 = 54.5 ± 5.9 kms-\ Kem = 306.1 ± 5.0 kms"^ and 
ifi = 0.007 ± 0.004 from the Rfi, 7 = 32.8 ± 9.7 kms"\ 
Kem = 312.6±6.4 kms"^ and ip = 0.001±0.007 from the H7 
line. The errors are obtained after we artificially increased 
the error-bars on the individual measurements such that the 
fit reduced was 1. The initial reduced '^^^ 1-6 for the 
H/3 line (45 d.o.f.) and 2.2 for the H7 one (33 d.o.f). The 
fit of the radial velocity curve of the two Balmer Hues pro- 
vides values consistent within la for ip and Kem, and within 
2a for 7. The latter is expected to be the same for all the 
lines, being the radial velocity of the system center of mass. 
ifi and Kem instead are not a priori the same for different 
emission lines, as they could originate from different parts 
of the irradiated companion star, or in fact regions in the 
binary. 

3.3 Averaged spectrum 

In order to measure the FWHM of the emission lines, 
we averaged the spectra in the frame of the companion 
star. We assumed a circular orbit, shifting the lines by 
V = —KainSin{2n(j) + ip), with Lp — 0. We did not adopt 
a priori the Kem obtained from the radial velocity curves 
of the Balmer lines, since that could be affected by, e.g., 
asymmetry or variations in the line profile from the indi- 
vidual spectra (as for the HQ:, |Bassa et al.|2009| ). Moreover, 
the Balmer lines could originate in a different area than the 
weaker Fell ones and thus have a different Kem. Instead, we 
proceeded as follows: 

- we assumed a range of possible values of Kem (200< 
Kd <400 kms^^ in steps of 10 kms"^). For each Kem value 
we applied the circular orbital shift and averaged the spec- 
tra between phase 0.3 and 0.7 (the range where the weaker 
emission lines are clearly detected). 

- we measured the FWHM of the emission lines in the 
averaged spectra, fitting a Gaussian function to each line 



(see Section 3.1 1. An orbital shift close to the real one will 
result in narrower lines in the averaged spectrum. 

- we found the Kem that results in the narrowest lines; for 
each line, we plotted the width measured from each average 



spectrum against the Kem of the corresponding orbital shift. 
The uncertainty on the width was large with respect to its 
variation among different values of Kem. Nonetheless, each 
line displays a trend indicating a minimum width at Kem ~ 
300 kms~^. The result is consistent with that from the radial 
velocity curves in Section [3. 2| 



The spectrum in Figure |4] is the average of individual expo- 
sures comprised between phase 0.35 and 0.7 (where all the 
emission lines are visible) for Kem ~ 300 kms~^. Together 
with II/3 and H7, the averaged spectrum highlights the weak 
Fe (or He) emission lines. 

No absorption lines appear in the averaged spectrum. As 
for the individual spectra, a cross-correlation of the average 
spectrum with the standard template spectra does not give 
a match. Note that we have also constructed an averaged 
spectrum in the phase interval 0.8 to 0.2 and 0.9 to 0.1, but 
we did not detect absorption lines from the non-irradiated 
face of the companion star. 

Table [1] presents the parameters of the best-fitting Gaussian 
to each line in the averaged spectrum. Each line was fitted 
for FWHM, normalisation and for the offset A — Aq of the line 
centroid A with respect to the rest-frame wavelength Ao. As 
we corrected the spectra for the orbital motion only, the lines 
in the average spectrum are still shifted with respect to their 
rest-frame wavelength by the systemic radial velocity 7. The 
wavelength offset of the H/3 and the H7 lines is consistent 
on a 2a level with the values of 7 derived from the radial 



velocity curves (see Section 3.2 1. The weighted average of 



the messures of 7 from the radial velocity curves and from 
the offsets of the Balmer lines is 7 =43.8±3.6 kms-\ The 
offsets of the Fell lines are not consistent with this value, 
although they agree with each other on a 3a level. In order to 
understand this difference, we have tested the dependence 
of the measured line offsets in the averaged spectrum on 
the choice of Kem and p in the orbital motion correction. If 
the Fe 11 lines originate in different regions of the companion 
star, in fact, the two sets of lines will be associated with a 
different Kem and ip. We found that the offset changes slowly 
with Kem, but is sensitive to variations in p>. The offset of 
the Fe 11 lines in the average spectrum is consistent with the 
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Figure 4. Spectrum of the counterpart to EXO 0748—676. This spectrum is an average of th irty-one 900 s-long spectra collected in 
the phase interval 0.3 and 0.7, corrected for the companion star orbital motion (see Section 3.3 i. The feature at '~5577 A is due to non 
perfect subtraction of the night sky OI line at that wavelength. The two prominent emission lines are from H7 (4340.465 A) and H/3 
(4861.327 A). Four weaker lines are marked in the spectrum; the position of the first two is consistent with Hel lines, at 4921.929 A and 
5015.675 A, but also with Fe II hues at 4923.92 A and 5018.44 A. The third hne is consistent with the Fen line at 5169.03 A, which 
favors an interpretation of the first three lines as an Fell triplet. The fourth weak line can be identified with Mgl at 5183.604 A. The 
insert shows a zoom-in of the region of the Fe II triplet and the Mg I line. 



a) 



Ao(A) 



FWHM (kms-i) Normalisation (A^) Ofl'sct (kms^i) 



H7(4340.465) 


267±10 


21.9±0.71 


33.1±8.08 


H/3(4861.327) 


282±7 


29.9±0.7 


42.5±6.9 


Fe 11(4923.92) 


246±28 


5.6±0.55 


3.5±12 


Fe 11(5018.44) 


366±50 


5.7±0.61 


-50±21 


Fe 11(5169.03)1 


410±36 


8.4±0.66 


20±15 


Fe 11(5169.03)2 


311±41 


6.9±0.77 


-15±11 


b) 










FWHM (kms-i) 


V sin jcm{kms" 


-1) Ofl'set(kms-i) 


Fell group 


304±20 


255±22 


-15±11 



Table 1. a) Results from the fitting of individual emission lines in the average spectrum of EXO 0748—676 (see Figure |4]l with a 
Gaussian function. The wavelength of each line in the rest frame Aq (column 1) is frozen in the fit. The last column in the Table indicates 
the velocity offset of the fitted line centroid respect to Ao. There is evidence for the presence o f a n arrow spike on the red wing of the 
Fe ll(5169.03) line: we fitted it without masking the spike Q-) and masking it {^) (see Section |3.3| for a discussion), b) Results of the 
combined multi-Gaussian fit of the three Fell lines, fitted together forcing the same offset and FWHM in kms"^. The spike affecting 
Fe 11(5169.03) is masked. The V sin iem is measured by artificially broadening Arc lines with a Grey profile until reaching the FWHM 
measured for the EXO 0748—676 lines. This accounts for the instrumental resolution profile. The line smearing due to orbital motion 
during the integration time is also taken into account. 
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Ao(A) 



FWHM (kms-l) Normalisation(A^) Offset (kms-l) 



Mg 1(5183.604) 172±39 



2.2±0.4 



8±16 



98±39 



Table 2. Best parameters from the Gaussian fit of the faint emission Une identifiable with Mgl (5183.604 A) in the average spectrum 
(see Figure |4|. The last column reports the intrinsic line width i)siniom, measured as for the Fell lines (see caption Tab.[l]l. 



measure of 7 from the Balmer lines if the orbital motion of 
their source region is shifted in phase by — 0.03< ip < —0.05. 
As we previously pointed out, the FWHM of the lines (both 
H7, H/3 and the Fell group) is not affected by changes in 
Kcm of a few tens of kms~^ around ~300 kms~^. Unlike 
the offset, it is also not sensitive to changes of a few percent 
in ip. In other words, our measure of the FWHM is not af- 
fected by a possible small displacement of the Fell source 
region with respect to the source region of the Balmer lines. 
All the lines in the average spectrum have the same FWHM 
at a 2a level, with the exception of the reddest Fe line, Fe 
(5169.03 A), which is significantly broader than Fe (4923.92 
A) at the >3(T level. Careful inspection shows a narrow faint 
emission line on the red end of the 5169.03 A line (Figure 
|4|. The presence of this second peak artificially broadens 
the Gaussian function used to fit the Fe (5169.03 A) line, 
as the fitting routine tries to account for both peaks with 
the same Gaussian. If this narrow peak is masked or fitted 
with an additional line, the width of Fe (5169.03 A) becomes 
consistent within la with the other lines. It must be noted 
that a second Gaussian introduced to fit the faint emission 
near the Fe (5169.03 A) line helps the fit but is not signifi- 
cant itself (the normalisation is approximately equal to the 
uncertainty on it). After we have verified that the Gaussian 
fit to the Fe 11 lines separately provides FWHMs and offsets 
in agreement with each other, we fitted the triplet together 
in order to reduce the uncertainty on the parameters. The 
three Fe lines were forced to have the same FWHM in km s~^ 
and the same velocity offset. The faint emission on the red- 
dest Fe line was masked to avoid artificial broadening. We 
obtained a FWHM=304 ± 22 kms~^ (reduced = 0.99, 
994 d.o.f. ). 



3.4 Measuring the intrinsic broadening of the 
emission lines 

The observed FWHM of the emission lines is determined 
by the intrinsic line width broadened by the instrumental 
resolution profile and by line smearing produced by the mo- 
tion of the companion during the integration time of one 
spectrum. Intrinsic saturation effects can also contribute to 
the width of the emission lines, significantly affecting strong 
lines like the Balmer ones ( Osterbrock||1989 1. On the other 
hand, if the weak Fen lines originate on the companion 
star, their intrinsic width is expected to be dominated by 
the rotational broadening usiniem ( Gray|19 92[). In order to 
measure the latter from the Fell triplet, we have artificially 
smeared and broadened the arc spectra with different val- 
ues of V sin Jem . The v sin iom resulting in smeared arc lines 
as broad as the FWHM of the Fe 11 lines is a measure of the 
rotational broadening. We proceeded as follows: the phase 
of each EXO 0748—676 spectrum participating in the av- 
erage in Figure [4] was ascribed to one arc spectrum. We 
smeared the arc spectra by 2-KTKcm cos{2-Kcj>) / P where T is 



the duration of one exposure (900 s) and P the orbital pe- 
riod. The smeared arcs were then averaged, each one with 
the same weight as that of the corresponding (same phase) 
EXO 0748—676 spectrum in the average of Figure |4] In this 
way we have simulated the effect of the smearing due to the 
orbital motion on the average spectrum. We then broadened 
the smeared average arc spectrum with a Grey profile for dif- 
ferent values of wsiniem, measuring each time the resulting 
FWHM in kms"^ We find that the measured FWHM of 
the Fell lines corresponds to i;siniom= 255± 20 kms~^. 



3.5 Evidence for Mgl in emission 

The average spectrum displays a weak emission line consis- 
tent with Mgl at 5183.604 A. The line was not detectable in 
the single or trailed spectra, but a fit with a single Gaussian 
function indicate that the line is significant on a 5(t level 
in the average spectrum. Two other Mgl lines, at 5167.321 
A and 5172.684 A, belong to the same multiplet with Mgl 
(5183.604 A). The position of Mg 1 (5172.684 A) is consis- 
tent with that of the spike on the red wing of Fell (5169.03 
A) (see Fig. |4]|, although the detection of the latter is not 
statistically significant when the prominent Fen line is in- 
cluded in the fit. Mgl (5167.321 A) could be present but 
it can not be resolved from Fell (5169.03 A). The FWHM 
of the Mgl (5183.604 A) line is 172±39 kms~\ which cor- 
responds to an intrinsic width of «siniem= 98± 39 kms~^ 



(measured as we did for the Fe 11 triplet in Section 3.4 1. This 



value is smaller than what we obtained from the Fe 11 triplet. 



3.6 Doppler tomography 

We employed emission line Doppler tomography to map the 
observed emission features. Here we used the modulation 
Doppler tomography method of Stccghs (2003), using the 
same code that was employed in the Bassa et al.^ (2009| 
study. The observed spectra were first phase binned using 
30 orbital bins. The underlying continuum was subtracted 
using a polynomial fit to line-free regions. In Figure [5] we 
show the resulting tomograms for the H/3 line and the two 
Fen lines at 4923.92 and 5169.03 A. The H7 line recon- 
structions were effectively identical to H/3 and we compare 
the two Fe 11 lines given their apparent difference in intrinsic 
width. All emission lines reconstructions are consistent with 
line emission from a region near the mass donor star with 
no evidence for any extended emission from a residual ac- 
cretion disc or stream. We find that the line flux is strongly 
modulated, as expected for an origin on the irradiated face 
of the companion, and in phase with the observed contin- 
uum modulation. Good fits are achieved close to a reduced 
of 1. 

The Doppler maps confirm the key attributes we derived in 
the previous sections using Gaussian fitting. We find that 
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Figure 5. Tomographic reconstructions for the H/3 and for the 
Fell lines at 4923.92 A and 5169.03 A. The top-left panel in 
each case shows the observed line profile as a function of the 
orbital phase while the top-right panel shows the reconstructed 
data from the converged maximum entropy solution. The bot- 
tom panels show the Doppler maps resolving the emission of each 
line in a Yx—^Y plane. The constant and the phase-dependent 
contribution to the flux are shown in the left and right quadrant 
respectively (see |Steeghs|2003) . The color scale indicates the frac- 
tional amplitude of the variable map in per cent. A Roche lobe 
model and ballistic stream trajectory is plotted for an assumed 
system mass ratio of q=0.25 and radial velocity semi-amplitude 
of the companion star K2 =410 kms~^. 



the radial velocity amplitude of the emission ranges from 
300-350 kms~^, consistent with the Kcm values from Sec- 
tion |3.3| as well as past estimates of K2 . We show in Fig- 
ure [5] the expected location of the donor Roche lobe and 
gas stream using a mass ratio o f q=0.25 (c onsistent with 

and the work 



2006 



the NS mass range favoured by Ozel 
of Muiioz-Darias et al.|2009 see also Bassa et al.|2009 \ and 
7^2=410 kms~ . The latter is the lower limit to K2 obtained 



by Bassa et al. ( 2009 1 which is a more stringent constraint 
compared to what we derived in the previous sections (see 
the discussion in section [5|. Again, the hne at 5169.03 A 
appears somewhat anomalous with a significantly more ex- 
tended emission distribution, reflecting its observed width. 
Furthermore, the center of the emission is shifted slightly 
off the Roche lobe, consistent with the apparent phase shift 
mentioned previously. 

In summary, our Doppler tomograms confirm that the nar- 
row emission lines are consistent with an origin on the irra- 
diated face of the mass donor star, with the exception of the 
emission feature at 5169.03 A which is both extended and 
shifted with respect to the other lines. This is in agreement 
with the broadening e ffect observed in the Gaussian fitting 
of 5169.03 A (Section 3.3 1. We will discuss this further in 
Section m 



4 SUMMARY OF THE RESULTS 

The results presented in the previous sections can be sum- 
marised as follows: 

- The g'-band hght curve of the source spans a range 
of 1 magnitude. The profile is variable, but the errors on 
the individual magnitudes are large: significant variability 
at a 3(7 level is detected within the phase interval of 0.007 
starting from phase 0.349, with an amplitude of 0.35±0.08 
magnitudes (Figure |2] top panel ). The morphology of the 
light curve in the g'-band is similar to that of the r'-band 
light curve reported by [Bassa et al. (20091, albeit the latter 
is not variable. 

- The spectra show emission lines from H/9 and H7 and 
three weaker lines consistent with Fell (4923.92 A) or He I 
(4921.929 A), Fell (5018.44 A) or Hel (5015.675 A) and Fell 
(5169.03 A). Since there is no He I line at the position of the 
last weak emission feature, we favour an Fe 11 group interpre- 
tation. A weak emission line consistent with Mgl (5183.604) 
is also detected at a 5a level in the average spectrum. 

- The orbital motion displayed by the emission lines is 
consistent with them originating close to the surface of the 
companion star facing the NS. The phasing of the contin- 
uum variation is in sync with the strong orbital modulation 
of the emission line strength, as expected if both are mainly 
produced by the irradiated companion. Evidence of an ir- 
radiated companion star was previously found, both during 



the X-ray outburst phase (Pearson et al. (20061, Muiioz- 



Darias et al. 20091 and at the beginning of the quiescent 
period ( [Bassa et al.|2009[ ). 

- The broad disc component at the base of the emission 
lines that was detected in the spectra from [Bassa et al.[ 
( 2009 1 is not detected here. There is no evidence for an 



accretion disc either from the line profiles or the Doppler 
tomography. 
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- The fitting of tiie radiai veiocity curves of tiie H7 and 
H/3 lines provides 1 a consistent measures of tlie radial ve- 
locity amplitude of the emission lines Kom, whose weighted 
average is 308.5±3.9 kms~^. The minimum width of the 
emission lines in the average spectrum is obtained with an 
orbital- motion correction with Kem ~300 kms~^, in agree- 
ment with the radial velocity curves estimate. This mea- 
surement provides a lower limit on the radial velocity semi- 
amplitude of the companion star K2, which is fully consis- 
tent with the results oflMuiioz-Darias et al. (20091. A more 



stringent lower limit on K2 is reported by Bassa et al. ( 2009 1 
(see discussion in section [5|. 

- The systemic radial velocity 7 was measured from the 
radial velocity curve traced by the H7 and H/3 lines and 
from the wavelength offset of the same lines in the orbital- 
motion subtracted averaged spectrum. The resulting values 
of 7 are in agreement within 2a, their weighted average being 
43.8±3.6 kms~^. The offset of the Fell triplet in the average 
spectrum, instead, is inconsistent with this value of 7. A 
phase shift of -0.03/-0.05 of the source region of the Fell 
lines with respect to the binary axis can account for this 
discrepancy (see Section [5|. The Doppler maps also show 
that the emission of the Fe 11 lines is shifted, particularly for 
the one at 5169.03 A. Our measure of the systemic radial 
velocity is in agreement with the work of |Mufioz-Darias et al.| 
( 20091), but is more than 5cr apart from the value obtained 



by 



Bassa et al. 



( [2009| ), 7 =78.6±3.9 kms"\ Even taking 
into account that the value of 7 is affecte d by the accuracy 
of the wavelength calibration (0.023 A in Bassa et al.||2009 



and 0.1 A for our data) the two measurements differ by 
more than 4(t. Similar mismatches using VLT data are not 



unprecedented (see Orosz et al.||2011 1 . 

- With the exception of the weak Mg I, the emission lines 
are very broad. The FWHM of the Fen lines indicates a 
V siniem= 255± 22 kms~^, if the broadening results from the 
rotation of the companion star alone. The Mgl (5183.604 
A) line instead provides a much lower 11 sin of 98± 39 
kms~^. 



5 DISCUSSION 

Using medium resolution VLT F0RS2 spectra of 
EXO 0748-676, ranging from 4222 to 5701 A, we 
found emission lines indicating that the companion star 
is still heated by irradiation, as was observed soon after 



due to a change in a pulsar wind (see below) or to vari- 
ability in the X-ray radiation, although the latter is hard 
to explain given the small variations observed in the X-rays 
( Degenaar et al.|2011 Dfaz Trigo et al.|2011 1. Furthermore, 
theoretically one does not expect the X-ray flux of the glow- 



ing cooling NS to vary on short timescales (Brown et al 



20021. A comparison with the R-band light curve presented 
Bassa et al. (12009) shows that the morphology of the light 



curves is similar, albeit the g'-band one is variable. 
Although we can not rule out a disc contribution to the 
optical continuum, the non-detection of disc line emission 
is unusual in comparison with typical observations of qui- 



escent LMXB (e.g. Marsh et al. 1994 Torres et al. 20021 



The disc-instability model that can explain the short out- 
burst of typical X-ray transients ( |Smak|[l971[ |Dubus et al.| 



2001 1 predicts that the disc is present even during quies- 



cence and that an outburst is triggered by a sudden rise of 
viscosity. The lack of disc lines in EXO 0748—676 suggests 
instead that the mass transfer has dramatically dropped and 
might have even stopped since the observations performed 



by Bassa et al. ( 2009 1 at the beginning of the quiescent 



phase. On the other hand, the disc-instability model might 
not apply to the case of EXO 0748—676 as it can not easily 
explain the 20-year long outburst that the source underwent. 
The trigger of tens of years long outbursts is not well under- 
stood and might have to do with variations in the envelope 
of the companion star, rather than with instabilities of an 
accretion disc ( [Remillard fc McClintock||2006 K 
In an attempt to obtain constraints on the NS mass from 
the emission lines, we determined the rotational velocity 
obtained from the FWHM of the Fell triplet observed in 
the spectra. If those lines are produced on the companion 
star surface, their width is expected to be dominated by 
the effect of rotational broadening v sin iem ( Gray|1992 1 . As 
also emerged from the tomographic reconstructions, the line 
emission does not arise from the full inner hemisphere of the 
companion star, but concentrates close to LI point. For this 
reason the v sin iem obtained from the Fe 11 lines would rep- 
resent a lower limit to the projected spin velocity of the star 
vsini. Due to the shift between the center of light of the 
Fe emission lines and the center of mass of the companion, 
the radial velocity semi-amplitude of the emission lines Kem 
also provides only a lower limit to K_d. The combination of 
the lower limits on Ko and v sin i provides in turn a lower 
limit on the NS mass. From the radial velocity curves and 
from our average spectrum we measured K2 > 308.5±3.9 



the onset of the quiescent phase (Bassa et al. 20091. We kms (see Sections 



3.2 



and 



3.3 1. From the Doppler maps 



detected H7, H/3 and a set of Fell lines, which is unusual 
but not unprecedented in LMXBs ( [Marsh et al]|1994| |. We 
also detect a faint Mg l line. The source of the irradiation is 
most likely the X-rays associated with the cooling compact 
object (Degenaar et al. 20091. The slow rate of the NS 
cooling that has been observed in EXO 0748—676 with 
XMM , Chandra and Swift (Degenaar et al. 2011 Dfaz 



(which are not influenced by asymmetries in the line profiles, 
see 
kms" 



Bassa et al.||2009| l we obtain a higher limit, K2 >355 



et al 



An even more stringent limit is provided in Bassa 
(2009|, K2 > 405 kms~\ This was obtained from 



Doppler tomography of the Ha line, which is expected to 
form closer to the center of mass of the companion star with 
respect to higher ionization potential lines such as the H/3 



Trigo et al. 20111 is in agreement with the fact that the 



companion star is still significantly irradiated, after more 
than one year of quiescence. 

The g'-band light curve of EXO 0748—676 shows evidence 
for variability between different orbits, in the phase interval 
of 0.007 starting from phase 0.349. The variability ampli- 
tude is 0.35±0.08 magnitudes. This implies a change in the 
fiux heating the companion star, which could potentially be 



and H7, providing a radial velocity closer to Kd (e.g. Har- 
laftis et al.(T999l|Unda-Sanzana et"aL][2006| . 



The constraints on the mass of the NS in EXO 0748-676 
that can be obtained combining the results of [Bassa et al.| 
( [2009 1 with our measure of v siniom show that, if the width 
of the Fe lines in our spectra is indeed setting a lower 
limit to usini, then the mass of the compact object in 
EXO 0748-676 should be >3.5 Mq. This is more than the 
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value of ~3 Mq which is considered to be an upper bound 
to the mass of a NS. A BH interpretation would be in con- 
flict with the detection of thermonuclear X-ray bursts from 
the source, which establish that the compact object is a NS 



(Gottwald et al. 19861. Given the fact that a bursting BH 



or a 3.5M0 NS would be against theories about those ob- 
jects, we have investigated possible scenarios to solve this 
incongruity. 

Our mass constraint is based on the assumption that the 
width of the Fe lines provides a measure of the projected 
rotational velocity of the companion star wsini, but that 
might not be the case. An asymmetric irradiation of the 
companion due, e.g., to the shadowing of a residual accre- 
tion disc, could induce currents in the atmosphere trying to 
re-distribute the heat ( |Kirbiyik| 19821 [Wolff et al.|2002| ). The 
velocity of the current would add to the width of the emis- 
sion lines. Nevertheless, this velocity has to be of the order 
of the total line broadening in order to give a significant 
contribution, i.e. a few hundreds kms~^. This is more than 



the sound speed on the surface of a late type star (Frank 
et al.|2002| |. 



Another possibility is that the Balmer lines, originating 
closer to the companion star surface, are broadened by sat- 



uration effects (Osterbrock 19891, while the Fen triplet. 



slightly phase shifted, is produced in a stream of residual ac- 
creting matter whose velocity causes the observed FWHM. 
This interpretation requires a bit of fine tuning of two broad- 
ening effects and it does not fit well with the lack of any 
evidence for an accretion disc. 

The large width of the Fe lines could also be due to a blend 
with unresolved lines of other atomic species. This is sug- 
gested by the detection of a faint (5cr significance) line con- 
sistent with Mgi at 5183.604 A. As mentioned in Section 
|3.5| the latter belongs to a Mgl multiplet including lines at 
5167.321 A and 5172.684 A, both very close to the position 
of the reddest detected Fell line. Mgl (5172.684 A) might 
be identified with the peak next to Feii (5169.03 A) that 
we have masked when measuring the FWHM of the Fell 
line. The peak is not significant in the average spectrum, 
but the anomalous extension and shift of Fen (5169.03 A) 
in the Doppler maps suggest that this narrow feature might 
be real and present in all the spectra, altering the look of 
the line in the Tomographic reconstruction. On the other 
hand, Mgl (5167.321 A) could be present, but it is not re- 
solved from Fell (5169.03 A). While a blend with Mgl can 
explain the width of Fell (5169.03 A), a blend with He I can 
broaden the other two lines of the Fell triplet. As in the 
"stream" scenario, the FWHM of H/3 and H7 would instead 
be determined by saturation effects. The downside of this 
sketch is that it requires some fine tuning of the width and 
intensity of the lines of different species in order to produce 
consistent FWHM and offsets among the three Fell lines. 
Other lines from Mgl and He I (at 4387.928 and 4471.68 A) 
that are within the wavelength coverage of our observation 
are not detected, so we have to meet the further requirement 
that only the lines that overlap with Fe H are strong enough 
to be observed. Although the Mgl (5183.604 A) detection 
supports this scenario, we can not confirm it or rule it out 
with the resolution and signal-to- noise of our observations. 

For this reason, we have considered another possible 
situation, where the lines are produced in a wind of mat- 
ter evaporating from the companion star, e.g. due to the 



action of a pulsar wind and/or the X-ray heating. The emis- 
sion lines, originating in the wind of evaporated material, 
are broadened due to the wind velocity. The Balmer lines, 
produced closer to the companion star surface, display the 
orbital motion of the irradiated secondary, while the Fell, 
rising from further out in the wind, display a slight phase de- 
lay in the orbital motion as the evaporated material starts to 
trail the companion star as it moves outwards due to angular 
momentum conservation. The evaporation of material could 
also explain the lack of evidence for emission lines from an 
accretion disk, as it would quench the accretion. Moreover, 
the presence of a pulsar wind in addition to the X-ray heat- 
ing might be responsible for the observed variability in the 
g'-band light curve, which is possible but rather difficult to 
produce in terms of variability of the X-ray fiux. At last, this 
"Black-widow like" scenario is consistent with the extended 
distribution of the line emission in the Doppler maps. The 
weak point of this outline is that it may be difficult to pro- 
duce the Mgl (5183.604 A), which is narrower than the other 
lines. The weakness of the line did not allow the production 
of Doppler maps to verify if the source region of this feature 
is narrow or extended and if it is the same as the Balmer 
or Fe 11 lines. Moreover, it has to be noticed that the pul- 
sar in quiescence has not been found in the radio waveband 
(private communication M. Burgay). 



6 CONCLUSIONS 

We performed phase-resolved medium-resolution optical 
spectroscopy with VLT/FORS2 on the eclipsing LMXB 
EXO 0748—676, one year and two months after it returned 
to quiescence. We found evidence for variability in the g'- 
band light curve and we found that the spectra only present 
emission lines coming from the region close to the compan- 
ion star, with no line contribution from a disc. This is un- 
usual compared with typical X-ray transients and could be 
related to the mechanism that triggers long duration out- 
bursts such as those displayed by the source, which is not 
yet understood. 

The g'-band variability can be explained by a varying X-ray 
fiux or by a variable pulsar wind. X-ray observations show 
little variations in the fiux, as expected from a cooling NS. 
The H/3, H7 and Fell emission lines in the average spectrum 
are remarkably broad. Their FWHM is not likely to be dom- 
inated by rotational broadening alone as this would lead to 
a compact object mass inconsistent with a NS. The width 
of the lines can be explained with bulk gas motion, e.g. the 
lines could be produced in a wind of material evaporated 
from the companion star envelope due to X-ray heating, 
possibly in combination with a pulsar wind. It is also possi- 
ble that some of the lines are blends of umesolved features. 
Despite extensive efforts to search for spectral features using 
multiple methods, we detect no absorption features from the 
mass donor star. Deeper observations with better resolution 
might allow to resolve blends and detect more emission lines 
that could provide a reliable lower limit to the companion 
star rotational broadening and, in turn, to the NS mass. 
Moreover, there is still the chance to detect absorption lines 
from the non-irradiated face of the companion star, which 
would finally provide a mass estimate for this elusive NS. 
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